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ABSTRACT 

The richness of dynamical behavior exhibited by the rotational states of various solar system objects 
has driven significant advances in the theoretical understanding of their evolutionary histories. An 
important factor that determines whether a given object is prone to exhibiting non-trivial rotational 
evolution is the extent to which such an object can maintain a permanent aspheroidal shape, mean¬ 
ing that exotic behavior is far more common among the small body populations of the solar system. 
Gravitationally bound binary objects constitute a substantial fraction of asteroidal and TNO popu¬ 
lations, comprising systems of triaxial satellites that orbit permanently deformed central bodies. In 
this work, we explore the rotational evolution of such systems with specific emphasis on quadrupole- 
quadrupole interactions, and show that for closely orbiting, highly deformed objects, both prograde 
and retrograde spin-spin resonances naturally arise. Subsequently, we derive capture probabilities for 
leading order commensurabilities and apply our results to the illustrative examples of (87) Sylvia and 
(216) Kleopatra asteroid systems. Cumulatively, our results suggest that spin-spin coupling may be 
consequential for highly elongated, tightly orbiting binary objects. 


1. INTRODUCTION 

The rotation of most natural satellites of the solar sys¬ 
tem is synchronous with their orbital periods (Peale 1999; 
Murray & Dermott 1999). Beginning with the seminal 
study of Darwin (1879, 1880), this fact has long been 
understood to be a result of dissipative planet-satellite 
interactions (MacDonald 1964; Kaula 1964). With con¬ 
siderable emphasis placed on describing the rotational 
evolution of the Moon (Goldreich 1966), the wide-spread 
applicability of spin-orbit synchronization as the end- 
state of tidal evolution was generally accepted as of half 
a century ago (Liu & O’Keefe 1965; Peale & Gold 1965). 
However, the discovery of Mercury’s a-synchronous spin 
(Pettengill & Dyce 1965; see also Margot et al. 2007) 
overturned the conventional understanding of tidal evo¬ 
lution of the time (Peale 1988). 

Shortly after the observational revelation, it had be¬ 
come clear that Mercury’s 59-day rotation period is 
a natural consequence of spin-orbit coupling, an effect 
that qualitatively alters the isolated process of tidal de¬ 
spinning (Golombo 1965; Goldreich & Peale 1966, 1968; 
Golombo & Shapiro 1966). In turn, the added insight 
has been instrumental to understanding the dynamical 
evolution of natural satellites (Wisdom et al. 1984; Wis¬ 
dom 1987, 2004; Marcialis & Greenberg 1987; Dermott 
et al. 1988; Black et al. 1995; Kouprianov & Shevchenko 
2005; Tiscareno et al. 2009; Cuk & Nesvorny 2010) as 
well as the dramatic history of the Moon itself (Mignard 
1979, 1980; Burns & Matthews 1986; Touma & Wisdom 
1994, 1998; Ganup 2004; Aharonson et al. 2012; Cuk & 
Stewart 2012; Makarov 2013). 

Spin-orbit coupling arises from gravitational torques 
exerted onto a permanently deformed satellite by the 
central object. Such modulation is particularly conse¬ 
quent when the torques accumulate coherently (i.e. when 
there exists a near-rational relationship between the or- 
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bital and rotational frequencies) and strongly depends on 
the degree of triaxiality of the satellite (see Ch. 5 of Mur¬ 
ray & Dermott 1999). Therefore, exotic rotational states 
are naturally expected to be more prevalent among the 
solar system’s small body populations, since typical ob¬ 
ject shapes become increasingly more irregular as mass 
is decreased^ (see for example Thomas 2010; Busch et al. 
2011 and the references therein). 

Evidence for complex spin dynamics among sub¬ 
planetary objects exists within the current observational 
census. Specifically, Hyperion (a small satellite of Sat¬ 
urn) as well as Nix and Hydra (minor satellites of Pluto) 
are observationally inferred to rotate chaotically (Klavet- 
ter 1989; Showalter & Hamilton 2015). Besides the 
highly unusual satellite figures, in the case of Hyperion, 
chaos is facilitated by an eccentric orbit (Wisdom et al. 
1984; Wisdom 1987), whereas in the Plutonian system, 
irregular motion stems from perturbations due to an¬ 
other satellite, namely Gharon (Gorreia et al. 2015; Gol¬ 
dreich & Peale 1967). While the families of spin-orbit 
interactions at play in the aforementioned examples are 
distinct, these systems share a common feature: the cen¬ 
tral objects in both cases are almost perfectly spherical. 

Not all satellites in the solar system orbit spherical 
bodies. In particular, the nearly three-decade old inquiry 
of Weidenschilling et al. 1989 - “Do asteroids have satel¬ 
lites?” has now been definitively and positively answered 
(see Merline et al. 2002 - “Asteroids do have satellites”). 
In fact, asteroidal and trans-Neptunian binaries comprise 
a non-negligible fraction of the overall small body popu¬ 
lation. Among near-Earth asteroids, ~ 15% are thought 
to be binaries, with a similar fraction corresponding to 
small {R < 10 km) main-belt asteroids (Bottke & Melosh 
1996; Pravec et al. 2006). Among larger asteroids (10 km 

^ As an example, consider the fact that the asphericity of Mer¬ 
cury {R = 2440 km) is of order {B — A)IC ~ 10“^, whereas that of 
Hyperion (R = 135 km) is (B — A)fC ~ 0.25 (Murray &; Dermott 
1999; Wisdom 1987). 
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< i? < 100 km), the binary fraction is probably some¬ 
what smaller (e.g. ^ few %), similar to that of the dy¬ 
namically hot component of the Kuiper belt (Merline et 
al. 2002; Richardson & Walsh 2006; Noll et al. 2008). 
The binary fraction among the cold classical population 
of the Kuiper belt is considerably more enhanced and is 
thought to be of order ~ 25% (Stephens & Noll 2006). 

There are a few different formation channels that may 
be associated with binary populations occupying differ¬ 
ent parts of the solar system. A process most relevant to 
the formation of near-Earth and small main-belt aster- 
oidal binaries appears to be rotational fission, facilitated 
by radiative torques i.e. the YORP effect (Bottke et al. 
2000, 2002; Scheeres 2007; Walsh et al. 2008; Pravec et al. 
2010). The primary mode of formation of larger binary 
asteroids is likely associated with (sub-)catastrophic im¬ 
pacts (Durda 1996; Durda et al. 2004; Doressoundiram 
et al. 1997; Michel et al. 2001; Scheeres et al. 2002). In 
trans-Neptunian space, dynamical capture (Goldreich et 
al. 2002; Astakhov et al. 2005; Noll et al. 2008), colli¬ 
sions (Pan & Sari 2005; Canup 2005; Stern et al. 2006; 
Schlichting & Sari 2009; Leinhardt et al. 2010), and fis¬ 
sion during gravitational collapse (Nesvorny et al. 2010) 
are believed to dominate binary formation. 

Irrespective of the exact formation mechanism, highly 
deformed objects smaller than R < 100 km can be found 
throughout the solar system. Thus, a considerable frac¬ 
tion of all minor object binaries may comprise triax- 
ial satellites that orbit permanently deformed primaries 
(Marchis et al. 2006). Real life examples of such systems 
include the triple asteroid systems (87) Sylvia and (216) 
Kleopatra (Ostro et al. 2000; Brown et al. 2001; Marchis 
et al. 2005; Descamps et al. 2011). These systems are 
comprised of large (Rp ~ 200 km) primary objects, or¬ 
bited by minuscule {Rg ~ few km) moons. The orbital 
separation of the moons is of order a ~ few Rp and the 
longest dimension of the primaries is roughly twice that 
of the shortest dimension. The tidal spin down timescale 
associated with these satellites is of order tens of thou¬ 
sands of years (Goldreich & Sari 2009). 

How is the rotational evolution of the small moons 
affected by the highly irregular shapes of the central 
objects? Qualitatively speaking, the gravitational po¬ 
tential of a strongly triaxial body harbors a substantial 
quadrupolar component. As a result, during the course of 
tidal spin-down or spin-up, the rotation rate of the moon 
may become commensurate with the evolution frequency 
of the quadrupolar part of the potential, which in turn 
varies with the spin of the primary. In other words, ro¬ 
tational evolution of multiple small-body systems can be 
subject to spin-spin coupling. 

In this work, we will explore the spin-spin coupling 
effect quantitatively. We note that observational char¬ 
acterization of multiple asteroid systems is a relatively 
young field (the first triple asteroid, (87) Sylvia was dis¬ 
covered a decade ago - see Marchis et al. 2005), and di¬ 
rect measurements of the rotation rates of minor moons 
do not yet exist. Accordingly, the primary goal of this 
paper is to provide broad predictions regarding yet un¬ 
characterized modes of rotational evolution, rather than 
to perform a detailed study. Correspondingly, simplicity 
is favored over realism throughout the performed calcu¬ 
lations. 

The paper is organized as follows. In section 2, we 
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calculations, a refers to the physical orientation of the satellite, 
while (p corresponds to the rotation of the primary. / denotes the 
true longitude of the satellite. In a non-inertial (primed) coordinate 
system correlating with the primary, a and / are replaced by the 
angles b and 9 respectively. 

analytically derive the governing equations for low-order 
spin-spin resonances, and consider their stability. In sec¬ 
tion 3, we compute the associated probabilities of cap¬ 
ture. In section 4, we present illustrative applications 
of the obtained results to the multiple asteroid systems 
mentioned above. We conclude and discuss our findings 
in section 5. 


2. RESONANT TORQUES 

The coupling function we aim to derive here is one 
describing the rigid body dynamics of an aspherical sec¬ 
ondary (satellite) of negligible mass, orbiting a rotating 
aspherical primary (planet/asteroid/KBO). For defini¬ 
tiveness, only the secondary’s spin is assumed to evolve 
under the influence of the primary’s gravitational po¬ 
tential. Furthermore, we assume that the spin vectors 
of the objects and the corresponding principal axes are 
normal to the orbital plane. In other words, the formu¬ 
lation of the problem is restricted and planar. More¬ 
over, motivated by the nearly circular and planar orbits 
of (87) Sylvia and (216) Kleopatra systems (Descamps 
et al. 2011; Fang et al. 2012), we shall assume null eccen¬ 
tricities. 


2.1. Geometrical Setup 

The geometrical setup of the problem is depicted in 
Figure (1). For a considerable fraction of the following 
derivation, we shall work in a rotating reference frame 
centered on the barycenter of the primary and 
synchronized with its spin, (j}, such that the a:'-axis is al¬ 
ways aligned with the long axis of the primary ellipsoid. 
The radial line joining the secondary’s barycenter posi¬ 
tion {xQ,yQ) and the origin has the length r and makes 
an angle 9 with the x'-axis. The angular position of the 
secondary in the rotating frame (x', y') can be translated 
to that in an inertial frame (x, y) by adding the primary 
orientation angle to 0. As such, their sum is simply in¬ 
terpreted as the secondary’s true longitude f = 9 -\- (j). 

Following Hut (1981), we treat the secondary as a 
sphere of radius Rg, and represent its quadrupole mo¬ 
ment with two diametrically opposed point masses, mq, 
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placed a distance Rs away from the barycenter. In this 
formulation, the difference in the principal moments of 
inertia is interpreted as: Bs — As = 2mqRl (Murray & 
Dermott 1999), keeping in mind that As and Cs corre¬ 
spond to the long axis and the axis or rotation of the 
secondary, respectively. 

The positions of the masses mq in the rotating frame 
are arbitrarily taken to be and defining 

the distance between them and the primary’s barycen¬ 
ter as ri and r 2 respectively. The angle made between 
the a:'-axis and the line connecting {x'i,y[) to {x'^^y'^) is 
labeled as 4’ while a = ■0 -|- (^ is taken to be its inertial 
counterpart. 


2 . 2 . Derivation of the Torques 

With all relevant angles defined, let us now calculate 
the torque exerted on the secondary by the primary. We 
begin by expressing 

(a;o, 2 /o) = (rcos 6 »,fsin 6 l), ( 1 ) 

and 

{x[,y[) = (f(ecos '0 -|- cos 9), f{e sin tjj + sin 0 )), ( 2 ) 

where e = i?s/f <C 1. Additionally, we have 

ri = \/r^ + Rl + 2fRs cos {4 — 9). (3) 


Correspondent expressions for (a;^, j/ 2 ) r^ can be ob¬ 
tained by replacing by V' ~ tt in the above expressions. 

The torque on the secondary due to rriq at [x'i,y'4) is 
given by 

Ti = -mq{x[ - x'Q,y[ - j/o,0) x (4) 

where lA is the gravitational potential of the primary. 
To second order in the quantity (Rp/r), the primary’s 
potential can be approximated by MacCullagh’s formula 
(Murray & Dermott 1999): 

Grup G(Bp Cp — 2Ap)x‘^ 


_l_ G{Cp + Ap — 2Bp)y'^ 

where Ap,Bp, and Cp are the primary’s principal mo¬ 
ments of inertia^. Substitution of equations (1), (2), (3) 
and (5) into equation (4) yields the following expression 
for the orbit-normal component of the torque (the other 
components being zero by construction): 


1 = ^-yyyX 

8^3 (^£2 _|_ 2 ^ _j_ 2 £ cos(0 — 4)) 

( - 42e(Mp - Bp) sin(26») 

— 6e (2e^ — l) {Ap — Bp) sin(2^) 

— 6 (e^ -I- 1 ) {Ap + Bp - 2Cp) sin{9 - 4) 
+ Srupf'^ {^A + 3e^ -I- l) sin(0 — 4) 

— 6e{Ap + Bp — 2Cp) sin(2(0 — 4)) 

+ IdnipAe (e^ -I- l) sin(2(0 — 4)) 


^ Higher order expressions have been worked out by Ashenberg 
(2007) and the references therein. 


-|- Snipr'^e'^ sin(3(0 — 4)) 

— 15(Mp — Bp) sin(30 — 4) 

— 3e^(Mp — Bp) sin(0 — 34) 

— 3 (l3e^ — 1 ) (Mp — Sp) sin(0-I-V')). (6) 

As already mentioned above, the equation for T 2 can sim¬ 
ply be obtained by replacing 4 in the above expression 
by — TT. 

While complete, the above expression for the torque is 
rather cumbersome. Consequently, rather than working 
with the full expression for the total quadrupole torque 
Tq = Ti -I- T 2 , it is sensible to expand the torque as a 
power series in the small parameter e and Fourier de¬ 
compose the results. This procedure allows for a more 
straightforward study of the individual resonant harmon¬ 
ics. We take this approach below, limiting ourselves only 
to zeroth and hrst order resonances. 

2.3. Zeroth Order Resonances 

Setting mqC^ = {Bs — As)/2A in the above expression, 
expanding to leading order in e and switching to angles 
in the inertial frame, we obtain: 

( 0 ) ^ SGmpjBs - As) sin( 2 (/ - a)) 

9 2f3 

5{Ap Bp — 2Cp) 

Ampf^ 

_ 2\G{Bp - Ap)iBs - As) sin( 2 (a - </>)) 

16f5 

lObGiBp - Ap){As - Bs) sin(4/ - 2{a + 4)) 

16f5 ■ 

The first harmonic in the above expression is responsible 
for conventional spin-orbit coupling (Goldreich & Peale 
1966). In fact, if we take the primary to be completely 
spherical {Ap = Bp = Cp), the above equation reduces 
to the well-known expression for torque exerted onto an 
aspherical satellite by a point-mass (Danby 1988). 

The other two harmonics in the expansion are respon¬ 
sible for spin-spin coupling. In particular, both terms 
give rise to 1:1 spin-spin resonances. Note, however, that 
while the second term librates when the rotation rate of 
the secondary is commensurate with that of the primary 
(i.e. d ~ (/)), the third harmonic requires the rotation 
rates to have opposite signs. In a regime where the or¬ 
bital frequency is well separated from the rotational fre¬ 
quencies (i.e. / 0), this harmonic resonates when the 

rotation rates of the primary and secondary are equal 
and opposite (i.e d ~ —0). In other words, in the limit 
where / can be considered to be a slowly varying an¬ 
gle, the second term in equation (7) governs the prograde 
1:1 spin-spin resonance while the third term governs the 
retrograde 1:1 spin-spin resonance. 

Incidentally, here a direct analogy between spin-spin 
resonances and mean-motion orbit-orbit resonances can 
be drawn. The critical angles of mean motion resonances 
also comprise combinations of differences between the 
quickly-varying mean longitudes and the slowly varying 
secular angles, namely the longitude of perihelion and 
the ascending node (see for example Ch .8 of Murray & 
Dermott 1999). Mean motion resonances are generally 
labeled according to the coefficients in front of the mean 












4 


longitudes in the critical arguments. Consequently, in 
this work, we follow the established convention in label¬ 
ing the spin-spin resonances. 

2.4. First Order Resonances 

Following the same procedure as for zeroth order 
torques, the Fourier decomposition of the first order 
torques take the form 

( 1 ) _ 21 e G{Bp - Ap){Bs - A) sin(Q - 2(j) + f) 

9 - 8^5 

, 21eG{Bp-Ap){Bs-As)sm{3a-2<j)-f) , , 

+-§-5-■ (8) 

Both of these terms govern spin-spin resonances. In par¬ 
ticular the first term corresponds to the 1:2 resonance 
while the second corresponds to a 3:2 resonance. Note 
that all resonant terms obey D’Almbert rules: writing 
any critical argument as {jia + j2<P + jif)^ the sum of 
the coefficients j is always zero. This is because the setup 
of the problem is invariant under rotation of the inertial 
frame. 

We can in principle continue the expansion and derive 
additional higher order torques. However, the strength 
of the resonances rapidly decreases. As a result, here we 
shall stop the expansion at first order, leaving further 
elaboration for future studies. 

2.5. Stability of Spin-Spin Resonances 

As already mentioned above, for a spherical primary 
and a (nearly) circular orbit, the only stable end-state 
of tidal evolution is a 1:1 spin-orbit resonance, such as 
that observed for the Moon and numerous other plane¬ 
tary satellites. The presented calculation shows however, 
that this is no longer the case if both, the primary and 
the secondary are significantly non-spherical. Still, for 
an object to become permanently trapped in a spin-spin 
resonance, such a resonance must be stable (note that 
this is a separate issue from capture probabilities, which 
we consider below). As Goldreich & Peale (1966) point 
out, the stability of a resonance is assured simply by re¬ 
quiring the restoring quadrupole torque to exceed the 
tidal torque i.e. > Tj when the critical argument of 
the resonance approaches tt or 0 (depending on whether 
the resonance in question is prograde or retrograde). 

In the framework of the constant time-lag (CTL) tidal 
model (Mignard 1979, 1980; Hut 1981), the relevant tidal 
torque reads 

GraiRf 

Tt = -3fc2„T— ^ —(d-n) (9) 

where k 2 ^ = 3/2(1 -|- /1)“^ ~ 10“^(i?s/1000km)^ is the 
secondary’s Love number {p, is the dimensionless mean 
rigidity), n = / is the orbital frequency, and t is the 
(small) time by which the tidal potential is assumed to 
lag the perturbing motion (Peale 1999). Within the CTL 
model, the oft-quoted dissipation quality factor is related 
to the time lag by = 2r|d —n| (Efroimsky & Lainey 
2007). 

We note that although the CTL tidal model has been 
utilized extensively in various astrophysical contexts (e.g. 
Hut 1981; Correia & Laskar 2009 and the references the- 
rien), it only serves as an adequate approximation to 
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Fig. 2.— Level curves of the Hamiltonian (23). The red curves 
denote librating trajectories while the blue curves denote circulat¬ 
ing trajectories. The separatrix is shown as the black curve. The 
width of the resonance as well as the location of the fixed points 
are labeled. Note that both expand in time as oc The stable 
fixed points are labeled with filled dots, while the unstable fixed 
point is shown with an empty dot. Librational and circulational 
phase-space areas are also labeled. 



more realistic rheological models within a narrow param¬ 
eter range. As an example, for a body obeying Maxwell 
rheology, the CTL description can be employed provided 
that the viscosity is smaller than (Efroimsky 2015) 


4tt GplRl 
57 \a-n\' 


( 10 ) 


Strictly speaking, this implies that the following results 
apply to objects with viscosities orders of magnitude be¬ 
low that of ice. However, in the spirit of simplicity, here 
we shall follow the (formally incorrect) convention of ce¬ 
lestial mechanics and employ equation (9), leaving more 
complicated analyses for future work. 

Bearing in mind the form that quadrupole torques take 
from equations (7) and (8), the stability of spin-spin res¬ 
onances is ensured when 


3k2^TmpR^{a — n) 


fePf{Bp-Ap){Bs-As) 


3 ^2, mp 

k) 

1 ^ 

2 gl-p 

2^ Qs 

\Rp) 

AgAp 


where f is the resonant torque coefficient (generally of or¬ 
der ~ few - e.g. 21/8 in the case of first order resonances), 
p is the order of resonance and X = {B — A) fmR^. 

As an example, let us consider a ~ 100km sec¬ 
ondary orbiting a ^ 1000km primary, naively setting 
Q ~ 100 (Coldreich & Soter 1966). For such a binary, 
first order resonances require (AgAp) ~ 10“® for stabil¬ 
ity. For a physically smaller binary (i?s ^ 10km sec¬ 
ondary orbiting a Rp ~ 100km primary), the required 
value changes to (AgAp) ^ 10“^, thanks to the Love 
number’s dependence on the secondary’s radius. The re¬ 
quired (AgAp) further decreases by a factor of e for zeroth 
order resonances. Note, however, that if the secondary 
in question is a rubble-pile, tidal evolution proceeds at 
an enhanced rate as the mean rigidity gets replaced with 
peS —t \/pJY where Y ~ 10“^ is the yield strain (Col¬ 
dreich & Sari 2009). Consequently, the critical (AgAp) 
must also be somewhat higher for rubble piles. 
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Collectively, the above arguments suggest that only 
small deviations from absolute sphericity of both bodies 
are required for low-order spin-spin resonances to sta¬ 
bilize the secondary’s rotation against tidal de-spinning. 
As already mentioned above however, this fact alone does 
not guarantee that an encounter with a spin-spin reso¬ 
nance will result in capture. As a result, we shall consider 
spin-spin resonance capture probabilities in the next sec¬ 
tion. 

3. CAPTURE PROBABILITIES 

In the vicinity of a given spin-spin resonance, charac¬ 
terized by a critical angle (jia + j 2 (t> + jsf), all other 
quadrupole torques will average to zero over many cir¬ 
culation cycles. Consequently, it is sensible to consider 
averaged equations of motion where only a single reso¬ 
nant torque is retained: 

, CePG{Bp-Ap){Bs-As)sm{jia + j2(l) + j3f) 

c-sQ: + 

GwlR^g 

+ 3fc2^r ———{a — n) = 0 (12) 



1:1 spin-orbit resonance 




Changing variables to v = {jia + j 2 (t> + j^f), the aver¬ 
aged equation of motion simplihes to that of a damped 
pendulum: 

i)-I-sin(r;)-I-/Jd-I -7 = 0. (13) 

In the above expression, a circular orbit has been implic¬ 
itly assumed yielding the following constants: 

= A^tPn^Xp{Rplrf{Bs - A)/C« 

/3 = 2>k2^Tii?e^{mp/ms){msRl/Cs) 

7 = j2^(n-0). (14) 

Immediately, the argument made above about the sta¬ 
bility of resonances can be understood as a criterion for 
existence of librating trajectories (and by extension, the 
existence of the separatrix). In other words, if jql > uA, 
V can not reverse sign as a result of quadrupole torques, 
implying continuous circulation. Obviously, without a 
librational island in phase-space, capture into resonance 
cannot take place (recall however, that this is a weak 
criterion). Note further, that w is the natural libration 
frequency of the resonant angle, v. 

Following Henrard (1993), the equation of motion 
(13) can be understood as arising from a pendulum-like 
Hamiltonian 

^ cos(r;), (15) 


where I = v — ^/j3 is the conjugated momentum, per¬ 
turbed by dissipative forces: 


V = 


m 

dl 


i 



(16) 


Note that unlike mean-motion resonances of celes¬ 
tial mechanics, Hamiltonian (15) does not posses the 
d’Alembert characteristic (dependence of on the ac¬ 
tion; Henrard & Lamaitre 1983). 

There exists a general procedure for incorporation of 
dissipative effects into the Hamiltonian framework which 


Fig. 3.— A stroboscopic surface of section, showing the phase- 
space portrait of the inertial angle o;, corresponding to parameters 
characteristic of Remus, the inner satellite of the (87) Sylvia sys¬ 
tem. The parameters are chosen in accord to the observed proper¬ 
ties of the system: Rp/a = 0.2, € = 0.01, {Bp—Ap)/{nipR^) = 0.19. 
For the purposes of this figure, the same degree of triaxiality as that 
inferred for Nix is adopted i.e. {Bs — As)/Cs = 0.63. The three 
panels of the figure show resonances associated with the three har¬ 
monics that appear in equation (7). 

we shall utilize here (see Henrard 1993; Galley 2013; 
Tsang et al. 2015 and the references therein). Let the 
auxiliary system of equations 


^ = Fi i=F2, (17) 

where Fi and F 2 are dissipative forces and {Q,P) are 
initial conditions, be satisfied by the solutions v{Q,P,t) 
and I{Q,P,t)- These solutions can be interpreted as a 
transformation from the original variables, (w,/) to {Q,P) 
via V = V and 1 = 1. In other words, the Hamiltonian 
flow can be envisioned as an initial condition to the dissi¬ 
pative solution. The new Hamiltonian, K, that describes 
the system is then related to the old Hamiltonian, R by 

/C(g, P, t) = umQ, P, PA)), (18) 

while the new canonical time, cr is defined by the Wron- 
skian 


dt f dv dl dv di\ ~ 

^ ^ “ dP]^j ^ 


(19) 


where {} signifies the Poisson bracket. 

In the context of the spin-spin resonance problem, 
equation (16) implies that Pi = 0 while P 2 = —/3/. Con¬ 
sequently, we are presented with the following relation¬ 
ships: 

v = Q i=Pe-A, (20) 

Accordingly, the new Hamiltonian takes the form 

JC{Q.PA) = \[P^-^' -w'cos(Q) (21) 
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The Poisson bracket (19) evaluates to {u, = e 

defining a = e^*//3. 

In the current form, the Hamiltonian (21) yields equa¬ 
tions of motion where the momentum P is multiplied 
by an explicitly time-dependent function. Such formu¬ 
lation is awkward, and it is preferable to transfer the 
time-dependence to the coefficients. We can rescale the 
Hamiltonian (21) by a time-dependent function x(t) such 
that 

JC' = xm da' = ^ ( 22 ) 

X(t) 

Choosing x{t) = we obtain: 

IC'iQ, P,t) = ^{P- c{t)f - bit) cos(Q). (23) 

where c(t) = ( 7 // 3 )e^‘ and 5(t) = (we^*)^. Accordingly, 
the new canonical time is related to the original time by 
da' = The Hamiltonian (23) describes a pendu¬ 

lum, whose stable fixed point shifting upwards in phase 
space as c(t) oc and expanding as 1 / 6 ( 1 ) oc The 
corresponding phase-space portrait is shown in Figure 
(2). Both, the upwards shift and the expansion of the 
phase-space area, occupied by librational trajectories are 
necessary requirements for adiabatic capture to occur. 

The process of adiabatic capture can be understood 
intuitively as follows. Consider a circulating trajectory 
that initially resides above the pendulum’s separatrix. 
Eventually, such a trajectory will encounter the separa¬ 
trix and either enter the librational phase-space or drop 
down to the circulational phase-space below the separa¬ 
trix. Because the phase of the encounter is essentially 
random, the probability of capture is simply determined 
by the relative rates at which the trajectory is “invaded” 
by librational and circulational phase-space. In other 
words, Pcapture = Aiib/Atot- The phase-space areas oc¬ 
cupied by librating and circulating trajectories are given 
by Aiib = Ibwe^* and Adrc = 27r(7//3)e^‘ — Swe^*, while 
Atot = 27r(7//3)e^‘ -|- Swe^*. Consequently, 

Pcapture = 2 ^1 -|- ■ (^4) 

This expression highlights the argument made by Gol- 
dreich & Peale (1966) stating that capture into resonance 
requires the tidal torques to be dependent on the rota¬ 
tion rate of the secondary (i.e. /3 / 0). Conveniently, for 

the tidal model of choice, x/P = j 2 {n — 4>). 

Recalling that the orbit is assumed to be circular, equa¬ 
tion (24) confirms that capture into the 1:1 spin-orbit 
resonance is certain because j 2 = 0. However, capture 
into spin-spin resonance is only certain if 

( 25 ) 

This leads us to conclude that capture into spin-spin res¬ 
onances is likely only when the secondary’s orbit is in 
close proximity to the primary i.e. f/Rp ~ few. How¬ 
ever, even such a conhguration requires high asphericity 
of both bodies. We shall examine capture into spin-spin 
resonances for real-life illustrative examples below. 



Fig. 4. — Capture probabilities into prograde and retrograde 1:1 
spin-spin resonances as a function of dimensionless radius. Best- 
fit parameters are used whenever available, and as before, Nix’s 
triaxiality parameter is adopted for the satellite. The obtained 
probabilities are of order ~ few %, and are similar to that obtained 
for Mercury within the framework of the same tidal model. 


4. APPLICATIONS 

Let us now examine the presented calculations within 
the context of observed asteroidal systems, (87) Sylvia 
and (216) Kleopatra. To begin with, we shall explore 
the rotational phase-space portrait of Remus, the inner 
satellite of (87) Sylvia. The physical size of Remus is 
i?s ~ 7 km, meaning that it is approximately hve orders 
of magnitude less massive than the primary. Therefore, 
we can safely neglect its back reaction on to the primary 
and take both / and (j) as linear functions of time. 

The 5.18-hour rotation period of the primary is roughly 
6.4 times shorter than the orbital period of Remus. Thus, 
a simple way to construct a stroboscopic surface of sec¬ 
tion is to record the dynamical state (a, d) once every 
5 orbital periods or Remus, which in turn corresponds 
to 32 rotation periods of Sylvia. We note that another 
simple alternative is to section on the angle 9 and plot 
the non-inertial phase-space {ip, ip). 

Since for the system at hand Rg/f ~ 0.01, we have 
chosen to only retain zeroth order (in e) terms in the 
Hamiltonian. Meanwhile, the potential of the primary 
is modeled as that of a triaxial ellipsoid with dimen¬ 
sions 384 X 262 X 232 km (Fang et al. 2012). Although 
the physical shape of Remus itself is observationally un¬ 
constrained, a highly irregular figure is entirely plausi¬ 
ble given that numerous other solar system bodies (e.g. 
Janus, Epimetheus, Hyperion, Nix, Hydra, etc) are in¬ 
ferred to have [Bg — Ag)/Cs > 0.1. For definitiveness, 
here we adopt the same degree of triaxiality for the 
secondary object as that recently derived for Nix i.e. 
[Bg — As)/Cg = 0.63 (see Showalter & Hamilton 2015; 
Correia et al. 2015). 

The computed surface of section is shown in Figure 
(3) . The top and bottom panels depict portions of phase 
space occupied by prograde and retrograde spin-spin res¬ 
onances respectively, whereas the middle panel shows the 
conventional 1:1 spin-orbit resonance. Compared with 
spin-orbit resonance, the widths of spin-spin resonances 
are diminished by a factor of 


yVs-s ^ j Bp - Ap ^ Rp 

Ws-o y mpf2 ^ r 


( 26 ) 
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For the considered case, this factor is approximately 
^ 0.1, as can be confirmed from the Figure. An analo¬ 
gous surface of section can be made for the (216) Kleopa¬ 
tra system, however given the overall similarity in the 
parameters of the two systems, it would yield the same 
qualitative features. 

The chances that a de-spinning satellite will become 
permanently captured into a spin-spin resonance can be 
easily evaluated via expression (24). Retaining the same 
degree of triaxiality as above, we have computed the 
probability of capture for prograde and retrograde spin- 
spin resonances as a function of the scaled radius for 
(87) Sylvia and (216) Kleopatra. The obtained curves 
are shown in Figure (4) 

As can be immediately seen, at the locations of the 
inner satellites, capture probabilities are generally quite 
low, ranging from ~ 3% to ^ 10%, with the higher prob¬ 
abilities corresponding to the retrograde resonance (due 
to its somewhat enhanced width). Incidentally, these es¬ 
timates are commensurate with the ^ 7% probability ob¬ 
tained for Mercury’s capture into its current spin state, 
within the context of the same tidal model (Goldreich & 
Peale 1966). 

In the case of Mercury, the effect of core-mantle fric¬ 
tion enhances the chances of capture dramatically (albeit 
into the wrong resonances) (Peale & Boss 1977; Correia 
& Laskar 2009). In the same vein, it is noteworthy that 
uncertainties persist in the understanding of tidal evolu¬ 
tion of rubble-pile asteroids (Efroimsky 2015), and it is 
not inconceivable that there exist physical processes that 
can alter our simple estimates substantially. Barring the 
importance of such unaccounted-for effects however, our 
calculations suggest that spin-spin coupling will only af¬ 
fect the most exotic binary/multiple objects within the 
sub-planetary population. 

5. DISCUSSION 

Understanding the spin states of binary objects in the 
solar system is critical to disentangling their complex 
temporal evolution. While gravitationally bound small 
bodies are subject to both external and dissipative forces, 
the specifics of the rotational state of the binary can dic¬ 
tate the regime of operation of the associated effects (e.g. 
Goldreich & Sari 2009; Cuk & Nesvorny 2010). Goncomi- 
tantly, the full scope of spin dynamics of small objects 
has not been exhaustively explored, and remains an area 
of active research (Naidu & Margot 2015; Correia et al. 
2015 and the references therein). 

In this work, we have considered the rotational dynam¬ 
ics of highly distorted small satellites that orbit triaxial 
central bodies, with an eye towards identifying qualita¬ 
tive deviations from the conventional picture of tidal de¬ 
spinning and spin-orbit coupling. To this end, we have 
shown that when both bodies are sufficiently aspherical, 
spin-spin resonances ensue, raising the possibility that 
previously uncharacterized rotational behavior may be 
observed within the small body population of the so¬ 
lar system. For systems where primary spin and orbital 
frequencies are well separated, and satellites are negligi¬ 
bly small compared to their host bodies, spin-spin cou¬ 
pling takes on a particularly simple form, manifesting in 
two additional harmonics that librate when the absolute 
value of satellite rotation becomes commensurate with 
that of the primary. 


For the purposes of this work, we have deliberately lim¬ 
ited the scope of our calculations to circular orbits and 
only considered planar rotational motion. These simpli¬ 
fications allowed us to obtain a handle on spin-spin cou¬ 
pling from purely analytical grounds. However, in doing 
so, we have removed dynamical features from our de¬ 
scription, that can exist within the framework of a more 
complete treatment of the problem. One such effect is 
YORP (e.g. Scheeres 2007; Cuk & Nesvorny 2010), as¬ 
sociated with radiative torques exerted onto the bodies. 
Another effect that is well known to be consequent is 
chaos. 

Irregular motion in Hamiltonian systems arises from 
overlap of neighboring resonances (Chirikov 1979). 
Given that leading-order asynchronous spin-orbit reso¬ 
nance widths scale as oc \Je(B — A)IC, and high degrees 
of asphericity are required for spin-spin coupling to op¬ 
erate at a noticeable level, it is likely that for systems of 
interest, even small values of orbital eccentricity will im¬ 
merse the vicinity of the 1:1 spin-orbit resonance into 
an extensive chaotic sea (Wisdom et al. 1984; Laskar 
1996). Emblematically, this can be inferred from Figure 
(3), where the separatrix of the 1:1 spin-orbit resonance 
extends over nominal frequencies associated with the 3:2 
and the 2:1 resonances (see Murray & Dermott 1999). In 
contrast, it is possible that spin-spin resonances depicted 
in the same Figure will not be affected as much, since 
their equilibria are well removed from leading-order asyn¬ 
chronous spin-orbit resonances. Instead, one can spec¬ 
ulate that accounting for finite, but nevertheless small 
eccentricities will simply act to engulf the separatrixes of 
spin-spin resonances into thin chaotic layers, in accord 
with the conventional modulated pendulum paradigm 
(see Gh. 4 of Morbidelli 2002). 

A related point follows regarding attitude stability. By 
now it is well known that satellites undergoing chaotic ro¬ 
tation also exhibit irregular obliquity dynamics (Wisdom 
et al. 1984; Showalter & Hamilton 2015). On the other 
hand, satellites locked into stable spin resonances can re¬ 
tain stable spin-axis evolutions. Moreover, tidal forces 
typically act to damp obliquities, meaning that the end- 
states of quasi-periodic rotational and tidal evolution are 
configurations where the rotational and orbital vectors 
are nearly aligned. Therefore, for highly triaxial bodies, 
our assumption of planar rotational evolution is probably 
grossly violated in the vicinity of orbit-synchronous rota¬ 
tion, but is justified for the characterization of spin-spin 
resonances (which is the primary aim of this work). 

Having applied conventional adiabatic capture theory 
(Yoder 1979; Henrard 1982) to the problem at hand, we 
have calculated the probabilities for enduring spin-spin 
locking in the (87) Sylvia and (216) Kleopatra systems. 
While the obtained estimates are not negligibly low, 
they suggest that spin-spin resonant capture is generally 
unlikely, even for the closer satellites of these systems. 
Provided the limitations of the GTL tidal model em¬ 
ployed here, these estimates deserve to be reexamined 
within the context of a more realistic rheology. Never¬ 
theless, interpreting our results at face-value, we can 
expect that the type of rotational dynamics considered 
here will only impact bodies that are exceedingly closely 
orbiting, and will be somewhat uncommon within the 
overall sub-planetary census. On the other hand, unlike 



planets, small bodies are exceptionally numerous in the 
solar system, so by observationally probing ever higher 
magnitudes and recovering tighter binary orbits (see 
e.g. Pravec & Harris 2007), it may be the case that the 
number of known bodies for which spin-spin coupling 
is consequential, will eventually become sizable. Ac¬ 
cordingly, we predict that observed light-curves of such 
objects will reveal satellite rotation that is synchronous 


with the spins of the central bodies. 
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